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Electronic States and Superconducting Transition Temperature based on the 
Tomonaga-Luttinger liquid in P^Ea^UrO^-d 
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An NQR experiment revealed superconductivity of Pr2Ba4Cu70i5_a (Pr247) to be realized 
on CuO double chain layers and suggests possibility of novel one-dimensional(lD) superconduc- 
tivity. To clarify the nature of the ID superconductivity, we calculate the band dispersions of 
Pr247 by using the generalized gradient approximation(GGA). It indicates that Fermi surface of 
CuO double chains is well described to the electronic structure of a quasi- ID system. Assuming 
the zigzag Hubbard chain model to be an effective model of the system, we derive tight binding 
parameters of the model from a fit to the result of GGA. Based on the Tomonaga-Luttinger 
liquid theory, we estimate transition temperature (T c ) of the quasi-lD zigzag Hubbard model 
from the calculated value of the Luttinger liquid parameter K p . The result of T c is consistent 
with that of experiments in Pr247 and it suggests that the mechanism of the superconductivity 
is well understood within the concept of the Tomonaga-Luttinger liquid. 



KEYWORDS: 



Introduction 



Pr 2 Ba 4 Cu70i 5 _ i 
Luttinger liquid 



first-principles calculations, CuO double chain, superconductivity, Tomonaga- 



The newly discovered superconductor P^Ba^uyOis-a 
(Pr247) consists of semiconducting single-chain and 
metallic CuO double-chain besides the Mott insulating 
Cu02 plane. 1-3 It shows the transition temperature (T c ) 
of superconductivity up to T c ~ 20K in a moderate 
oxygen defect concentration range of 6 = 0.2 ~ 0.6, 
where 8 can be varied by controlling oxygens in the 
CuO single chains, and electrons are doped to CuO dou- 
ble chain. It is found that the main carrier is elec- 
tron from the measurement of the Hall coefficient at low 
temperature(T < 120K). An NQR (nuclear quadrupole 
resonance) experiment has confirmed the superconduc- 
tivity in the CuO double chain layers. 4 A spin-gap-likc 
behavior was also observed in the recent NMR experi- 
ment, 5 although it is not clear yet whether the spin gap 
is enough large. These experiments suggest that Pr247 
shows the possibility of novel one-dimensional(lD) su- 
perconductivity. 

In the previous papers, we address the ID double chain 
models to clarify the superconductivity of Pr247 from the 
theoretical point of view. 6,7 Based on the Tomonaga- 
Luttinger liquid theory, we analyze ID superconducting 
mechanism proposed generally by Fabrizio, 8 in which the 
spin gap plays an important role. In fact, the spin gap of 
the double chain model is found to become up to ~ 100K 
for a typical case and the ID superconductivity seems to 
be realistic in Pr247. 7 On the other hand, it is shown 
that the spin gap is very small at some parameter re- 
gion by using the density matrix renormalization group 
(DMRG). 9 Then, Nakano et al. proposed a spin fluc- 
tuation mechanism of the superconductivity for Pr247 
by applying the fluctuation exchange (FLEX) method 
to a quasi (ID) extended Hubbard model. 10 Berg et al. 
claimed that "CiS 3 / 2 " phase, in which the total spin 
mode is gapless but the half of the relative spin mode 
is gaped, is crucial to the superconductivity by using 



the bosonization theory. 11 However, these theoretical 
works employed the result of energy band calculation 12 
for not Pr247 but YBa 2 Cu 4 8 (Y124) in which CuO dou- 
ble chains are included with the same lattice structure, 
because the band calculation results of Pr247 has not 
been available. Although energy dispersions correspond- 
ing to the CuO double chain for the both compounds are 
expected to be almost equivalent, it is very difficult to 
avoid ambiguity of the model to be solved. 

In this article, we calculate the band structure of Pr247 
by the density functional calculation with the generalized 
gradient approximation of Perdew, Burke, and Ernzerhof 
by using the WIEN2k package. 13 We set the ID zigzag 
chain Hubbard model to be an effective model of CuO 
double-chain and obtain the tight-binding parameters of 
the model by fitting the energy dispersion. Using these 
parameters, we estimate T c of the model based on the 
Tomonaga-Luttinger liquid theory. It would give an ac- 
curate analysis of the electronic state of Pr247 and shed 
light on the mechanism of the novel quasi- ID supercon- 
ductivity. 

2. Band Calculation 

The calculations is done using the WIEN code, 13 
which uses a full potential-linearized augmented plane 
wave (FLAPW) calculation. The exchange correlation 
potential is treated by the generalized gradient approxi- 
mation of Perdew et al. The energy threshold to sep- 
arate core and valence states is -6 Ry. For the num- 
ber of plane waves the criterion Rmt (muffin tin radius) 
x -fi'max (for the plane waves) = 7. Self-consistency of 
the calculation is yield by using a set of 216 k points 
in the irreducible wedge of the Brillouin zone and the 
spin-orbit coupling is also considered in our calculations. 
The crystal structure is referred by the experiment val- 
ues obtained by Yamada et al. 2 To check the validity 
of our calculation, we also obtain the result of GGA for 
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Fig. 1. A band structure of Pr247 in k z = from non-spin polar- 
ized calculation. The Y point is at the corners, the E point in 
center, the X point at the midpoint of the horizontal edge, and 
the Y point at the midpoint of the vertical edge, 



Prl23(data not shown) and confirm that it agrees with 
previous study well. 15 

Figure 1 shows the band structure of Pr247, where di- 
ameter of circles represents weight of contribution from 
d x 2_ y 2 orbitals of Cu in the double chain, where the dou- 
ble chain is set to be along the 6-axis in our calculation. 
We find that seven bands cross the Fermi energy Ep on 
X — r line among them, two bands which have significant 
dispersion on T — Y line, correspond to the CuO double 
chain. It surely shows that CuO double chains are metal- 
lic and the dispersion along the c-axis perpendicular to 
Cu02 layer is very small as well as that of Prl23. Almost 
flat bands near and above Ep consists mainly of Pr's 4/ 
orbital. 

Fermi-surfaces of Pr247 are shown in Fig. 2 on the k x - 
k y plane with k z — 0. Six vertical lines in both sides 
present the Fermi-surfaces of single-chain and double- 
chain, where that of single-chain is denoted by the dashed 
lines, and the solid lines are corresponding to that of 
double-chain. Four closed Fermi surfaces are correspond- 
ing to CuC-2 plane. The hybridization between double- 
chain and Cu02 plane is very small, because the Fermi 
surfaces of the double-chain arc almost straight, though 
that of the single chain is warped. The result suggests 
that the band structure of the double chain is well de- 
scribed as a quasi-lD system. Although the existence 
of the Fermi surfaces of Cu02 plane seems to indicate 
Cu02 plane to be metallic, it becomes insulator because 
the electronic conduction in the plane is suppressed due 
to the so called Fehrenbacher-Rice state formed by the 
strong hybridization between Pr4/ and 02p orbitals. 16 
We note that it is not able to take into account this ef- 
fect within the GGA calculation. Therefore, we should 
not consider the position of the Fermi level obtained by 
GGA too serious. 

Next, we consider quasi- ID double chain Hubbard 
model as a minimum model reflecting the band struc- 



ture of Pr247. It contains intra-chain hopping integrals 
between Cu sites and inter-chain hoppings as shown in 
Fig. 3 . The Hamiltonian is given by 

H = tl Yl ( C i,m,a C i+hm,a + h.C.) 

.m,<7 Ci +' 2 , rn ,< 7 "r" a - c -) 

+ *3 X) ( C !m,<7 C i+3,m,<7 + h.C.) 
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+ U Y ( C I,m, CT C H-4,m,cr + h.C.) + 
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,m,a C i+3,m+l,a + h.C.) + (c i m (T Ci_3. m +i j(T + h 

i,m.i7 

where cj m a stands for the creation operator of an elec- 
tron with spin a at site i on m— th chain and n itm ^ = 
c i,m,a c i,m,a- Here, t„=i,2,3,4 is the hopping energy be- 
tween intra chain sites and ^=1,2,3,4 is that between 
the nearest-neighbor inter-chain sites as shown in Fig. 3. 
Here, U is the on-site Coulomb interaction parameter. 

In a noninteracting case (U = 0), the Hamiltonian eq. 
(1) yields the band energies as functions of wave numbers 
of k x and k y , 

E(k x ,k y ) = s a + 2ti cos(fc x ) + 2t 2 cos(2fc x ) + 2t 3 cos(3fc x ) 

+ 2i4 cos(4:k x ) + 2V\ cos(k y ) + AV2 cos(fc x ) cos(k y ) 
+ 4V 3 cos(2fca;) cos(ky) + 4V4 cos(3fc x ) cos(k y ). (2) 

We estimate the tight-binding parameters t n and V n 
of the model eq. (2) so as to fit E(k x , k y ) to the energy 
dispersion obtained by the GGA. The obtained values of 




Fig. 2. A fermi surface in k z = for Pr247 from non-spin polarized 
calculation. The V point is at the corners, the S point in center, 
the X point at the midpoint of the horizontal edge, and the Y 
point at the midpoint of the vertical edge, 




i 


U 






1 


-0.1227 eV 


-0.0203 eV 


2 


-0.4908 eV 


0.0060 eV 


3 


0.0876 eV 


-0.0097 eV 


4 


-0.0652 eV 


0.0026 eV 



Table I. The values of intra-chain transfer tj, and inter-chain 
transfer Vi for i-th neighbor sight. 



parameters are shown in Table I and the fitted dispersion 
is shown in Fig. 4 . The figure indicates that the obtained 
dispersion is in good agreement with the GGA result and 
our quasi-lD model is expected to simulate the electronic 
state of Pr247 very well. We also performed same calcu- 
lation for Prl24(data not shown) and found that £i/t2 
of Prl24 is smaller than that of Pr247. 17 



Fig. 4. The dispersion of the tight-binding model with the param- 
eters determined so as to fit the calculated band structure of 
Pr24, where the closed circles present GGA result. 

In the double chain model with two different Fermi 
points, low-energy excitations are given by a single gap- 
less charge mode with a gapped spin mode. 8,25-29 In 
this case, the SC and CDW correlations of the charge 
mode decay as ~ r~^p and ~ cos[2(kp 2 — k Fl )r]r~ 2K p , 
respectively, while the SDW correlation decays expo- 
nentially. Hence, the SC correlation is dominant for 
K p > 0.5, while the CDW correlation is dominant for 
K p < 0.5. Since the role of inter-chain couplings is con- 
sidered to stabilize the SC state at finite temperature, the 
electronic state in chain part is crucial to the creation of 
the SC. If K p > 0.5, the SC correlation is dominant in 
the chain and the total system becomes superconductor. 



3. Tomonaga-Luttinger liquid approach 

Next, assuming the spin gap to be large enough to 
exceed the transition temperature (T c ) of the supercon- 
ductivity, 7, 18 we consider T c by using the formulation of 
the Tomonaga-Luttinger liquid in cooperation with the 
mean-field analysis for quasi-lD systems. 20 It contains 
the Luttinger liquid parameter K p of chain part and cou- 
pling parameter A between chains, as shown below 

r \ ~i 1/(2-27) 

Tc = ^|-WVW^r , (3) 

2ira \ 4tt V 2 ' 2 J 2 J W 

where A = f^(2V? + 4F 2 2 + 4K, 2 + AVi), B{x,y) = 

T^eTy) ' an< ^ T = JR~' Here, vp is the Fermi veloc- 
ity. We approximate vp ~ vp*, where vp* — J" 1 ^ 2 
is an effective Fermi velocity constructed by two Fermi 
velocities; vp 1 and vp 2 in the double chain model. 18 ' 19 
Further, a is a short- wavelength cutoff parameter. 20,21 
This value related to the spatial extension of the quasi- 
particle, however, it is difficult to determine it without 
ambiguity. We should chose a £ 2, since the dimer state 
on the nearest neighbor sites is considered to composite 
main part of the quasiparticle. 22-24 



0.9 0.75 




Fig. 5. T c as a function of K p at n =0.5, 0.75 and 1.0 for 
a=2(dotted line), 3 (dashed line) and 4(solid line), respectively. 

In Fig. 5, we show T c as a function of K p for severl 
values of electron number per site n by substituting the 
fitting parameters of Pr247 into eq.(3). At first glance, 
the n— dependence of T c seems to be small, but we stress 
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that n— dependence of K p is not small. For example, near 
n = 1, K p becomes very small due to the umklapp effect. 
The n— dependence of K p for similar models has been 
already discussed in several works. 6 ' 7 ' 18 They indicate 
that K p has a maximum value almost around n ~ 0.7. It 
suggests that T c of Pr247 has also a peak near the same 
point. 

We roughly estimate K p of our model by the Hartree- 
Fock approximation; K p = / 1 U - , 18 ' 19 where this 

V ™F* 

approximation is valid only in the weak coupling limit 
and neglects the umklapp process. When U = 2eV and 
n = 0.75, we obtain K p ~ 0.79 and then T c ~ 12K 
for a— 2 and T c ~ 23K for a=3. These results are rea- 
sonable and consistent with the experiment results. 2 It 
also suggests that the novel mechanism of the quasi- ID 
superconductivity is realized in Pr247. 

4. Summary and Discussion 

In summary, we calculate the band dispersions of 
Pr247 by using the GGA. It shows that the electronic 
state of Pr247 is well described by the quasi-lD system. 
We also obtain the tight binding parameters by fitting 
the band dispersions to the quasi-lD zigzag Hubbard 
chain model. It would give an accurate analysis of the 
electronic state and a concrete basis to address the super- 
conductivity of Pr247. On the assumption that the spin 
gap is large enough to exceed T c , we estimate T c of the 
superconductivity in Pr247 by using the formula of T c for 
quasi-lD superconductor. The result is consistent with 
that of experiments and suggests that the mechanism 
of the superconductivity in Pr247 may be understood 
within the concept of the Tomonaga-Luttinger liquid. 

Finally, we would point out the possibility of the quasi- 
lD superconductivity in Y124, since the double chains 
are also contained in Y124. We estimate T c of the quasi- 
lD superconductivity in Y124 as ~ 40K for U = 2eV and 
a=2 by adopting the tight-binding parameters of Y124 
obtained by ref.10, where inter-chain hoppings is larger 
than those in Pr247 resulting higher T c . 

At this stage, the experimental evidence of quasi-lD 
superconductivity has not been explicitly observed, since 
the quasi-two-dimensional(2D) superconductivity(T" c ~ 
80K) in Cu02 plane masks the quasi-lD superconductiv- 
ity even if it exists. If we can suppress only the quasi-2D 
superconductivity, we expect the quasi-lD superconduc- 
tivity to be found in Y124 near 40K. In fact, the experi- 
ment of Zn-doped Y124 shows that T c have a plateau at 
T c ~ 25K in the case of doping rate being larger than 



2%. 30 If doped impurities destroy the superconductiv- 
ity of CuG-2 planes completely, we can expect that the 
plateau may be corresponding to the indication of quasi- 
lD superconductivity. Further both of theoretical and 
experimental studies in this regard is highly desired. 
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